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Abstract 

We report on the single crystal growth of the ferromagnetic Kondo lattice sys- 
tem CeRuPO using a Sn flux method. Magnetic susceptibility and electrical re- 
sistivity measurements indicate strong anisotropy of this structurally layered com- 
pound. They evidence that the magnetic moments order ferromagnetically along 
the c-direction of the tetragonal unit cell, whereas the crystal electric field (CEF) 
anisotropy favors the a6-plane. Therefore, CeRuPO presents the unusual case within 
rare earth systems, where the anisotropy of the interionic exchange interaction over- 
comes the single ion anisotropy due to the CEF interaction. 
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1 Introduction 



Intermetallic Kondo lattice systems liave attracted considerable attention in 
the last decades. While many Ce-based Kondo lattices show antiferromagnetic 
(AFM) ground states, only very few systems are known with ferromagnetic 
(FM) order and pronounced Kondo effects. Therefore, new FM Kondo lat- 
tice systems have to be found to study the underlying physics in more detail. 
Recently, we reported on the physical properties of polycrystalline CeRuPO, 
PP which seems to be one of the rare example of a FM Kondo lattice system 
(Tk ~ 10 K, Tq = 15 K). This compound, at the border between intermetallic 
and oxide compounds, crystallizes with the tetragonal ZrCuSiAs type struc- 
ture, which consists of alternating layers of RUP4 and 0Ce4 tetrahedra. [2] The 
results of LDA calculations indicate a strong two-dimensional (2D) character 
of the Fermi surface, therefore the anisotropy of the magnetic and electronic 
properties are of high interest. However, the synthesis of this material is diffi- 
cult, due to the high vapor pressure of phosphorus and the control of oxygen. 

The preparation of high quality single crystals is one of the key issues to 
study in detail the anisotropy and to understand the underlying physical in- 
teractions. In this contribution we will present the growth of CeRuPO single 
crystals from Sn flux and the study of the magnetic anisotropy by means of 
magnetic susceptibility x{T) and electrical resistivity p{T) measurements. 

Two different kinds of anisotropics determines the magnetic properties of a 
solid. The most obvious one is that related to the different x, y and z compo- 
nents of the magnetic moment. It induced e.g. an anisotropy with respect to 
the direction of an applied external field and is therefore comparably easy to 
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address. The second one is connected with different exchanges along different 
directions within the crystal structure, and can generally only be investigated 
with Q (momentum) dependent methods, like e.g. neutron scattering. In the 
present paper we shall only address the former one. In magnetic rare earth 
systems, this anisotropy is usually determined by the single ion anisotropy in- 
duced by the crystal electric field (CEF), which is generally quite strong and 
lead to a pronounced anisotropy of the magnetic susceptibility. In the case of 
Ce^^ in a tetragonal environment, the CEF splits the spin orbit J = \ mul- 
tiplet into 3 Kramers doublets, which can present either a strong Ising type, 
an easy plane or a Heisenberg behavior depending on the CEF parameters. In 
contrast, the contribution to this anisotropy originating from difference in the 
exchange parameters for the different x, and z component of the magnetic 
moment, is usually negligible in these rare earth magnetic systems. Our results 
demonstrate that in CeRuPO this contribution is larger than the effect of the 
CEF, and determines the orientation of the spontaneous FM moment. 

2 Experimental procedure 

Recently, we described the preparation of polycrystalline (PC) samples us- 
ing a Sn-flux method in evacuated quartz tubes. [1],[3] With quartz tubes 
the temperature range is restricted to 1000°C. It turned out that within this 
temperature it is not possible to get single crystals (SC) large enough for 
physical characterization. One reason might be that the liquidus temperature 
of CeRuPO, at the formerly used 82at% Sn-concentration, is above 1000°C. 
Therefore, we slightly changed the preparation route to achieve higher tem- 
peratures. First, we prereact the phosphorus with half of the Sn flux at 600°C 
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for 3 hours in an AI2O3 crucible, sealed inside an evacuated quartz tube. After 
quenching to room temperature, we opened the quartz tube and filled on top 
of the prereacted Sn-P the Ce, Ru02 and the other half of the Sn with an total 
molar ratio of 9:6:4.5:80.5 (Ce:P:Ru02:Sn). The AI2O3 crucible was put inside 
a Ta-crucible, which was closed under an inert atmosphere. This crucible was 
heated with a rate of 300° C per hour to 1500° C and kept at this temperature 
for 1 hour, subsequently the melt was cooled down to 900° C by moving the 
crucible out of the furnace with an average cooling rate of 4°C per hour. After 
the reaction the excess Sn was dissolved in dilute HCl, leaving platelets with 
a surface ~ 3 mm^. A typical SC platelet of CeRuPO is shown in Fig. [TJ The 
direction perpendicular to the surface corresponds to the crystallographic c- 
direction of the tetragonal unit cell, as shown by a Laue picture (see inset of 
Fig. [3]). As a side product, we also got SC of CeRu2P2 and Ce3Ru4Sni3, the 
physical properties of these new Ce-based systems will be discussed elsewhere. 

m 

Several SC were investigated with electron-microprobe and energy dispersive 
X-ray (EDX) analysis, no foreign phases could be detected. A typical EDX 
spectrum is shown in Fig. [2l All lines could be assigned to Ce, Ru, P, or O and 
a quantitative analysis reveals a stoichiometric Ce:Ru:P content. The presence 
of oxygen is confirmed by the Kq, peak at 0.53 keV; however, a quantitative 
analysis of this oxygen peak failed because the contribution to the spectrum 
is too weak, compared to the heavier elements. 

X-ray diffraction (XRD) patterns were recorded on a Stoe diffractometer in 
both transmission and reflection mode using a monochromated Cu-Kq, radi- 
ation (A = 1.5406 A). Magnetic measurements were performed in the tem- 
perature range from 2 — 300 K in a commercial Quantum Design (QD) mag- 
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netic property measurement system (MPMS) equipped with an RSO option. 
AC-transport resistivity measurements were carried out in a standard four- 
probe geometry using a commercial QD physical property measurement sys- 
tem (PPMS). 

3 Results and Discussion 

To determine the crystal lattice parameters of the SC, we performed XRD 
measurements in transmission mode on 4 powdered SC. The intensity of the 
resulting XRD pattern was low, due to the small amount of powdered SC. The 
lattice parameters refined by simple least square fitting, a = 4.027(3) A and 
c = 8.26(1) A are in good agreement with the reported single crystal data [2] 
and also with our PC samples. [T] However, the accuracy of this fit was poor, 
indicated by the larger errors of the calculated lattice parameters. Therefore, 
we performed a reflection experiment on a SC to determine the c parameter 
more accurately. The idea behind this experiment was to check the oxygen 
occupation between the Ce layers, because it was shown that the introduction 
of hydrogen in the closely related CeFeSi structure leads to a large change 
of the c parameter. [5] Fig. [3] shows the 00/ reflections of a CeRuPO single 
crystal measured in reflection mode, the bars (red) denote the theoretical peak 
positions from Zimmer et al, [2j the OOZ peaks are marked with a star. It is 
obvious that even at higher angles the observed 00/ peak positions coincide 
with the data from Ref. [2j. The calculated c from the higher order peaks is 
csc = 8.259(1) A, and is therefore only slightly larger than the values from the 
PC samples, where we found cpc = 8.256(2) A. Compared to the large differ- 
ence (~ 15%) of the c parameter going from clean CeFeSi to hydrogenated 
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CeFeSiH, we conclude no large difference in the oxygen occupancy between PC 
samples and SC in CeRuPO. In PC samples we have shown a stoichiometric 
oxygen occupancy using carrier gas-hot extraction, thus our present analysis 
evidences also a stoichiometric oxygen occupancy in the SC of CeRuPO. 

To study the magnetic anisotropy of CeRuPO, we have performed dc-magneti- 
zation measurements on a SC for ± c and if || c. A large anisotropy is 
evident from Fig. HJ where we show the magnetization measurements in the 
ordered phase at 2 K for both field directions. For H c the magnetization 
increases linearly with a slope of 1.1 ;Ub/T at small fields and saturates above 
fioHci = IT with an absolute value of fi'^^^ = 1.2 fi^- For H \\ c we observe a 
spontaneous magnetization at small fields, followed by a small linear increase 
(which might be due to small misalignment), saturating also above i^ci with 
a magnetic moment of /ig^^ = 0.43 ;Ub. In the inset of Fig. |l]we have plotted 
the well defined hysteresis loop, together with the virgin curve for H \\ c 
with a spontaneous magnetic moment of 0.3 fiB and a coercive field of He = 
0.1 T. For i7 ± c we could not observe any hysteresis. This behavior strongly 
suggests that CeRuPO is a collinear FM with the moments aligned along the 
c-direction, but with the easy magnetization direction perpendicular to the 
c-axis. The observed behavior is most likely due to a competition between the 
CEF anisotropy and the anisotropy of the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) exchange interaction with respect to the x, y and z component of 
the magnetic moment. A similar behavior was also observed and studied in 
detail by means of neutron and magnetization measurements on the FM Kondo 
system YbNiSn by Bonville et al. [6j On the one hand the magnetic exchange 
interaction (RKKY), which is responsible for the FM order, is larger for the 
z component of the moment than for the x and y component, yielding a 
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spontaneous FM ordering of the moments along the c-direction. On the other 
hand, the CEF anisotropy leads to a ground state CEF doublet with a larger 
saturation moment in the basal plane than along the c-axis. This scenario 
explains both the absence of hysteresis as well as the linear increase of the 
magnetization for if _L c, which is due to the rotation of the FM moments 
towards the afe-plane. Above Hd the FM moments are completely aligned 
along the afe-plane with the full moment of the CEF ground state doublet 
for this direction. An alternative scenario invoking a canted AFM ordering 
with AFM moments slightly canted away from the a6-plane can be excluded, 
because of the presence of inversion symmetry in the structure, which prohibits 
a macroscopic net FM canted moment. 

Now we discuss the magnetization data on the basis of the CEF model, whose 
Hamiltonian for the J = 5/2 multiplet of the Ce^"*" ion located in tetragonal 
point symmetry can be written as 



where 5^ are the Stevens coefficients and are the Stevens equivalent 
operators. [7\ The resulting three doublet eigenstates can be expressed with a 
single mixing coefficient rj |8] 



It is now possible to determine the ground state wave function by comparing 
the theoretical values of the saturated moment with the experimental data 



-f^CEF = B2O2 + -B4O4 + -B4O4 , 



(1) 




(2) 



(3) 



(4) 
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for the two different directions. We found that the ground state is the 
wave function and the excited states are the two Ti eigenstates, because for 
this configuration the saturation moment along z is fi^^^ = I^'l/^b ~ 0.43 /xb 
and /iga^t = fs'L/^B ~ 1.29 /iB along x, {ql = | is the Lande factor for Ce^"^) 
which is in very good agreement with the experimental obtained values for 
the two different field directions. If one assumes one of the as the ground 
state, the maximal value of n^^^ = ^^/Eghl^B ~ 0.96 /iB for i] = l/\/2 is 
well below the experimental /i^f^ and rule out such a configuration. From the 
temperature dependence of the entropy we have shown that the two excited 
Kramer's doublets are 6 and 30meV above the ground state. [1] Therefore, 
we can draw a preliminary CEF level scheme, shown in the inset of Fig. HI 
However, we can not determine the exact values of the Stevens coefficients 
from Eq. [H which would need a more sophisticated analysis and/or further 
experiments e.g. inelastic neutron measurements. 

The complex anisotropic behavior of CeRuPO is also reflected in the tem- 
perature dependence of the susceptibility x(T) measured along both crystal- 
lographic directions shown in Fig. [5l x(^) for if || c presents the typical 
behavior of a FM system. The ordering below Tc = 14K is clearly visible in 
a strong increase of the susceptibility, reaching a constant value below 5K. 
This increase is much weaker for fioH = 1 T. For H c the magnetic field de- 
pendence of x{T) is not so pronounced. At fiQH = 0.1 T a small peak appears 
below Tc, similar to what we have observed in PC samples. [1] The precise 
physical origin of this antiferromagnetic signature is still unclear, but we think 
that it is due to the complex anisotropy, because a similar feature was also 
observed in YbNiSn for a direction perpendicular to the FM axis (see Fig. 
7 in Kasaya et al. ^). We also have performed field cooled experiments for 
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both directions of the magnetic field at fioH = 0.1 T (sohd fines in Fig. [5]), but 
only for fioH ± c we could observe a small difference at low field, most prob- 
ably due to domain reorientation. At high temperature (T > 150 K) x~^(T) 
follows a Curie- Weiss law with 6^ ~ — 250K, 6^ ~ +4.2 K, and effective 
magnetic moments /icfr = 2.4 /xb for both directions, close to that of the free 
ion moment (2.54 /iB/Ce^'*', see inset of Fig. E]). The anisotropy of Gw confirms 
the conclusions drawn above, that the CEF anisotropy favors the spins to be 
perpendicular to the c-axis. The large difference between 9^ and 6^ implies 
a large positive B2 CEF coefficient, which is expected to result in a Fg CEF 
ground state, in agreement with our deduction based on the magnetization 
data at 2K. At lower temperature x^^(T) changes slope and the data can be 
fitted with a positive Weiss temperature of the order of Tq for both directions. 

The temperature dependence of the resistivity for the current fiowing perpen- 
dicular to the c-axis is plotted in Fig. El Measurements with current along 
the c-axis could not yet be performed because of the small thickness of the 
platelets. The overall temperature dependence is very similar to what we have 
reported on PC samples. [T] Three different temperature regions can be iden- 
tified: (I) Above 50 K the sample shows a linear behavior typical of a con- 
ventional metal. (II) Below 50 K there is a pronounced decrease, which is a 
distinct feature of a Kondo lattice system and can be explained by the onset 
of coherent Kondo scattering, due to the hybridization between the localized 
4/ and the conduction electrons. (Ill) In contrast to the PC samples, in the 
ordered state p(T) follows a power law dependence with p oc T^, indicating the 
presence of a Landau Fermi liquid. In the inset of Fig. [6] we have plotted this 
dependence as p vs. and the linear fit (red curve) gives p(T^) = po + AT'^, 
with Po = 5 pQcm and A = 0.11 fiQcraK^^. 
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The anomaly in the resistivity at the FM ordering temperature is visible in 
Fig. [7l where the temperature dependence of the derivative of the resistivity, 
dp/dT vs. T is plotted for two different samples with the magnetic field par- 
allel and perpendicular to the c-direction. We have shown previously on the 
PC samples that dp/dT{T, H) behaves similar to the temperature and field 
dependence of the specific heat and can be therefore taken as a direct measure 
of the development of the FM order. In the upper panel of Fig. [7] the mag- 
netic field is applied along the c-direction, whereas current is flowing within 
the a6-plane. For fiQH = 0, dp/dT present a sharp peak at T = 13K which 
broadens and shifts to higher temperatures with increasing magnetic field. For 
fioH = 7T the peak is situated at 16.5 K. The field dependence is stronger for 
field applied within the afe-plane (see lower panel of Fig. [7]). There the peak 
shifts from 13 K at zero field to 20 K for fiQH = 7T. An increase of Tq with 
magnetic field is expected for a FM, since the external field is supporting the 
exchange field. The stronger increase for field along the basal plane compared 
to field along the c-axis can be related to the higher saturation moment in the 
basal plane, which leads to a larger free energy gain in the external magnetic 
field. 

We now discuss the small difference in Tq of the SC compared to our previ- 
ous study on PC samples {T^^ = 15 K). [I] The Tc of the SC (T^c = 14 K) 
was determined from a Curie Weiss Fit of the x^^{T) data below 30 K. The 
anomaly in dp/dT{T) confirms this T^*-', if one takes into account the asym- 
metric shape of the peak (see Fig. [7]). However, the peak is broader in the 
SC data compared to the PC measurement, and therefore it is more difficult 
to determine Tq^ . The specific heat data for a SC gives similar results (not 
shown), the peak at Tq" = 14K is rather broad. The reduced ordering tem- 
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perature and the broader anomalies at Tq^ indicate that the crystallinity of 
the PC material is better than in the here presented SC. This is also evident 
from a smaller residual resistivity ratio for the SC Psook/Po*" ~ 30 compared 
to P3o^k/Po'^ — 50 in the PC sample. However, we can exclude a large differ- 
ence in the oxygen occupancy for SC and PC, as wc have shown above. Thus 
the broader transition at a lower Tc is likely due to a larger defect density, 
hke e.g. a P-0 substitution. 



4 Conclusions 



In summary, we succeeded in growing single crystals large enough to determine 
the anisotropy of CeRuPO by means of magnetic susceptibility and resistivity 
measurements. We showed that the spontaneous FM moments point towards 
the c-axis, while the saturation moment is larger in the basal plane. The 
latter property indicates an easy plane single ion CEF ground state, while 
the former property implies a stronger exchange for the z component of the 
moment than for the x and y component, this anisotropy overcoming that of 
the CEF (at low T). A more detailed analysis gives for the CEF ground state 
the Tq eigenstate, composed of | ± |) wave functions. In addition, we have 
presented detailed EDX and XRD results which show that the difference in 
the oxygen occupancy is small for PC samples compared to the SC. However, 
the SC have a slightly smaller ordering temperature {Tq^ = 14 K) compared 
to the PC samples (Tq*^ = 15 K). Further studies are necessary to determine 
the precise origin of this difference. 
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Fig. 1. (Color online) Photograph of a CeRuPO single crystal (length scale 1mm). 
The direction perpendicular to the sm'face is parallel to the crystallographic 
c-direction. 
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2. (Color online) Energy dispersive X-ray spectrum of a CeRuPO single crystal. 
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Fig. 3. (Color online) X-ray diffraction pattern of a CeRuPO single crystal, in re- 
flection mode. The bars denote the theoretical peak position [2j, stars indicate the 
00/ peaks. The inset shows a Laue picture of a CeRuPO single crystal, the surface 
is perpendicular to the c-axis of the tetragonal unit cell. 
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Fig. 4. Isothermal magnetization as a function of applied magnetic field at T = 2 K 
for H J- c (black symbols) and H \\ c (red symbols). The lower inset shows the 
well defined hysteresis curve for || c. The upper inset presents the proposed CEF 
energy level scheme. 
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Fig. 5. Magnetic susceptibility vs. temperature at fioH = 0.1 T (closed symbols) 
and IT (open symbols) for H \\ c (red) and H J- c (black). The solid line (blue) 
presents a field cooled experiment for _L c at 0.1 T. x^^{T) for both directions is 
shown in the inset, together with the Curie Weiss fits for T > 150 K. 
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Fig. 6. (Color online) Temperature dependence of the electrical resistivity of 
CeRuPO with the current perpendicular to the c direction. The low temperature 
part is shown in the inset as p vs. T^, below 5K p oc indicates the formation of 
a Fermi-liquid ground state. 
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Fig. 7. (Color online) Derivative of the resistivity of CeRuPO with respect to tem- 
perature plotted as dp/dT vs. T. for two different samples with -ff || c (a) and H J- c 
(b). The anomaly due to the onset of FM order is clearly visible as a sharp peak. 
It shifts to higher temperatures with applied magnetic field. 
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